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SECTION I 


FORTRAN PROGRAM LISTING 


The following pages show the FORTRAN program listings, Each 
program or subprogram is begun on a new page, except that the fluid 
properties package, Subroutine UDG1 through Function UDG9, has been 
listed as a single unit. 
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SECTION II 


tens a 


EXAMPLES OF USE OF PROGRAM 


1. DISCUSSION OF EXAMPLE A 


The example shown is the design of a singic-stage compressor and 
is based upon a design described in Reference 1, The input data cards 
are shown in the next subsection, and selected output data are shown in the 
following subsection. An analytic meanline rotor blade and an arbitrary 
meanline stator are used, so that with the aerodynamic calculations all 
three program sections are used. 


The original design was made using 15 streamsurfaces to describe 
the flow, but in the interests of eccnomy, only 8 have been used here. 
The same computing station geometry is used here as was used in the 
original design, and alternate streamsurfaces, as derived previously, 
are input to the analytic meanline section for the rotor geometry determi- 
nation, Streamsurfaces are specified from one station upstream of the 
leading edge to one station downstream of the trailing edge. Blade section 
specifications were also reproduced directly, and these specify a poly- 
nominal camberline blade with an uncam.bered leading edge. Output of 
NASTRAN data cards is specified by setting INAST = 4. Thus four-point 
averaging will be used, and no pressure load cards will be punched. 
Input to the interface routine includes a shape correction factor of 0.7, 
and a single deviation fraction curve, as was used originally. One degree 
of additional deviation at the trailing edge (over that cetermined by the 
basic deviation rule) is specified, and the point of maximum camber is 
given as 0.75 of the blade chord. The loss form indicator NLOSS is set 
to 4 for each station within the rotor blade, so that (in the aerodynamic 
analysis) loss coefficients will be determined as a function of the trailing 
edge values. Accordingiy, NLOSS is set to 1 for the trailing edge, anda 
loss coefficient distributior. is given. The results show the blade inlet, 
outlet, and camber angles as positive; they are set positive in this routine 
so that they may be used in conjunction with the deviation rule data. The 
additional deviation remains positive, but the total deviation is shown 
negative, as it is used in the aerodynamic section in conjunction with 
(generally) negative rotor blade angles. 


The input printout from the aerodynamic section follows, and it is 
seen that blade geometry for stations 4 through 9 are to come internally 
frorn the analytic meanline section. These data are listed with the rest, 
and show, for stations 5 through 9, the deviation "rising" from about -1.5 
degrees to the trailing edge values of about -13.5 degrees. No loss 
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coefficients are shown for stations 5 through 8; these will be generated 
internally from values at station 9, and the loss fraction curve shown as 

the final input data. A simple linear increase from leading edge to trailing 
edge is indicated. For the stator, whirl velocity is specified and NOUT2 

is set to 1 from one station upstream of the bladerow to one station down- 
stream. Thus data will be output for the arbitrary meanline blade section. 
For each of the stations within the stator, an isentropic efficiency varia- 
tion is given, relative to conditions at the rotor inlet. At the trailing edge, 
a relative total pressure loss coefficient is given, because for both rotor 
and stator a loss re-estimation is specified. As NEVAL is negative for 
stations 9 1nd 15, the new losses will not be incorporated into the iteration, 
but will ve computed after the solution is converged. One family of loss 
parameter/ diffusion factor curves is given in the data, and is called out 

for both blades. For the rotor, NDEL is set to 2, and the specified Prandtl- 
Meyer expansion angle varies from zero at the hub to -5 degrees (that is, 5 
degrees of compression) at the tip. As leading edge geometry is given, 

the incidence angles will be added to these values. NDEL is set to zero 

for the stator, so that no shock loss will be computed. Pressure load cards 
for the NASTRAN program are specified for both rotor ana stator, rour- 
point averages are specified for the rotor, three-point for the stator. 


The output from aerodynemic analysis shows details of the flow 
through the machine, and is very similar to that shown in Reference 1. 
The loss coefficient re-estimation shows a relatively high rotor loss 
coefficient, due principally to a less optimistic assumption regarding 
precompression of the inlet flow. 


The stator blade design follows, producing a blade that corresponds 
to the air angle distribution cetermined through it in the aerodynamic 
section. A uniform 3 degrees of incidence and 1 degree of extra deviation 
are specified. INAST is set to -3, specifying the output of NASTRAN data 
cards for the stator, without pressure load cards (which are specified for 
the aerodynamic section). Three-point averages are to be used. Details 
of the resulting blade section for the innermost streamsurface are shown. 


2. INPUT DATA CARD LISTING FOR EXAMPLE A 


The input data cards are shown on the following pages, 50 cards 
per page. No completely blank cards have been used. 
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3. SELECTIONS FROM OUTPUT PATA FOR EXAMPLE A 
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4. DISCUSSION OF }.XAMPLE B 


This example uses only the aerodynamic section of the program, and 
demonstrates the turbulent mixing calculation and the blade wake transfer 
calculation. The machine for which the computation is made is the so-called 
NACA Five-Stage Compressor that is described in References 2 and 3. A 
conventional multi-stage analysis calculation is set up, with a cumputing 
station at each blade edge (that is, there are two stations between each blade 
row). For each of the ten blade rows in the machine, fixed relative outlet 
flow angle and (two-dimensional) loss coefficient distributions are specified, 
The distrihutions are ones that describe well the performance of the blade 
rows at an operating point in the viscinity of the machine's design point, 
except that the l' ss coefficients do not reflect any variation in blade section 
performance with blade section location. In practical terms, this means 
that no additional rotor tip :osses are incorporated. The wake transfer analysis 
was developed to account for migration of boundary layer material and is 
therefore applicable to the diffusion loss component of the relative total 
pressure loss coefficients. However, the subject compressor operates 
transonically in all stages at design speed, and so the rotor loss coefficients 
in the input data generally include a shock loss component in addition to the 
diffusion loss component, The input data item that controls the wake transfer 
rate, CONTR, is set to 1.0. The value is probably lower than it should be, 
and this compensates to some extent for the fact that (for the rotors) the 
whole loss coefficient is being manipulated instead of just the diffusion loss 
portion thereof. The input eddy viscosity, CONMX, is set to 2.0 ft@/sec 
whicn is higher than the probable correct val:e, 


A listing of the input data cards is shown in the next sub section, and 
the corresponding output from the program follows, The program has 
computed the results with difficulty, although just a few more iterations are 
apparently required fer total convergence. With a final IVFAIL of 5 and final 
IFFAIL of 0, the result may be taken to be converged for most practical 
purposes. The following points should be noted when examining the results. 
When wake transfer calculations are made, the loss coefficients shown in 
the results differ from those in the input. When turbulent mixing calculations 
are made, the !oss coefficients shown in the results differ from those in 
the input and efficiencies computed for streamlines may take on apparently 
impossible values because both of these parameters are computed from 
conditions on ''streamlines",. 


It is of some interest to examine the effects of the wake transfer and 
mixing calculations upon the computed compressor performance, and shown 
below are some performance parameters for the four possible calculation 
modes involving (or not involving) wake transfer and mixing, The final 
station (outlet) total temperature difference between casing and hub stream - 
lines is shown as ina‘cative of the effects of the various calculations 


a ee ee ee ee ee ee ee 


upon the distributions of properties of the flow. All computations were made 
with identical data, except that the input data items NTRANS and NMIX were 
varied as shown. 


Total Pressure Isentropic (T case - Th‘ 
NTRANS NMIX Ratio Efficiency at outlet, deg. 
0 0 5.188 . 8313 40.9 
0 l 5.264 . 8341 3.9 
1 0 4.524 . 8073 259.4 
1 1 4.512 . 8086 194.8 


It is apparent that the wake transfer calculation can drastically effect the 
machine overall performance and the distribution of properties of the flow. 
Also, the mixing calculation had considerably less effect on the overall 
performance, but did modify the distributions appreciably. It should be 
recalled that the calculations were somewhat simplified in that the relative 
outlet angles and base (two-dimensionc!) loss coefficients were fixed 
throughout. 


Some final notes are given, as follows. 


1) Correlation of empirical and computed results is needed to 
establish appropriate vales for the wake transfer rate 
constant, CONTR, and the eddy inviscosity, CONMX. 


2) In order to correctly apply the wake transfer calculation to 
transonic blade rows the computer program shoulc be 
modified so that shock and diffusion loss coefficient components 
may be separately manipulated, 


3) If itis desired to use the wake transfer calculation when 
computing stations are located within blade rows, the program 
will need modification. Probably the only change would be 
to replace I-1 with I + NLI (I) (separately computed) wherever 
it occurs in the call statement starting on card $AR§. 225. 


4) Both the wake transfer and mixing calculations substantially 
increase computation time. For a fixed number of iterations, 
each of these options taken separately increases compvtation 
time by a factor of approximately three, and taken together by 
a factor of approximately six, In the case of the wake transfer 
calculation option, taken above or with the nixing calculation 
option, this is further compounded by the fact that it increases 
the number of iterations required for solution by a factor of 
two or three. 


5. 


INPUT DATA CARD LISTING FOR EXAMPLE B 
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